Lapse of transmission phase and electron molecules in quantum dots 
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The puzzling behavior of the transition phase through a quantum dot can be understood in a 
natural way via formation of the electron molecule in the quantum dot. In this case, the resonance 
tunneling takes place through the quasistationary (doorway) state, which emerges when the number 
of electrons occupying the dot reaches a certain "critical" value, N cr . Our estimation of this quantity 
agrees with the experimental data. The dependence of N cr on the dot's size is predicted as well. 
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One of the challenging problems in mesoscopic physics 
is the puzzling behavior of the transmission phase 
through a quantum dot, embedded in an Aharonov-Bohm 
ring. It was found in scries of experiments performed 
by the Weizmann groupie that all transmission am- 
plitudes through different resonant levels of a quantum 
dot are in phase. This necessarily implies an unexpected 
lapse in the evolution of the transmission phase between 
different resonant levels. In addition, it was found in re- 
cent measurements^ that this phenomenon takes place 
when the number of electrons inside the dot reaches a 
certain "critical" value (N cr > 15)2.. In spite of many 
publications addressed to these experiments no fully sat- 
isfactory understanding has been found yet^£. 

In this Rapid Communication we demonstrate that the 
observed phase-lapse behavior of the transmission am- 
plitude can be naturally explained by implying the for- 
mation of electron (Wigner) molecules inside quantum 
dots, proposed in recent publications^^. Moreover, 
this framework allows us to estimate N cr and then to 
determine how it is varied with a size of the dot. In 
order to explain our model in a proper way, wc first elab- 
orate the physical nature of the transmission phase in 
the case of noninteracting and interacting electrons. In 
particular, we concentrate on the role of the Pauli prin- 
ciple that prevents different conductance resonances to 
carry essentially the same internal wave function. This 
point represents a formidable obstacle for resolving the 
puzzling behavior of the transmission phase for different 
models of the quantum dot. We demonstrate, however, 
that this difficulty can be overcame in the context of the 
Wigner-molecule when an unstable state is developed in 
the middle of the dot. 

Let us consider the resonant tunneling through a quan- 
tum dot, represented by a potential Un{x), Fig. [T] The 
bottom of this potential can be moved by the plunger 
electrode, so that one observes the current sweeping 
trough different resonant states (E\) of the dot. We 
would treat this problem in the framework of a tunnel 
Hamiltonian approach. This approach is more trans- 
parent for evaluation of the transmission phase than the 
standard scattering theory, in particular, when the Pauli 
principle and the electron-electron interaction are taken 
into account. We introduce therefore the following tun- 
neling Hamiltonian: H = Hi + Hr + Hp + Ht, where 
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FIG. 1: (Color online) Resonant tunneling trough a quantum 
dot. j«i(H) are Fermi energies in the left (right) reservoir. The 
dotted lines show the potential Ud(x), needed for evaluation 
of the bound state wave functions in the Bardeen formula. 
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Here, aj_(oi,r) is the creation (annihilation) operator of 

an electron in the reservoirs and d k (dk) is the same op- 
erator for an electron inside the dot. For simplicity, we 
consider electrons as spin-less fermions. The term He 
denotes the Coulomb interaction between electrons in 
the dot and fi| [fir ] is the coupling between the states 
Ei(E r ) and E k of the reservoir and the dot, respectively. 
In the absence of magnetic field, all couplings fi are real. 

All parameters of the tunneling Hamiltonian ([I]) are re- 
lated to the initial microscopic description of the system 
in the configuration space. For instance, the coupling 
fi;(A is given by the Bardeen formula^ 

h 2 f <-» 

where 0/s(x) and Xl(r){ x ) are the electron wave functions 
inside the dot and the reservoir, respectively, and S;( r ) 
is a surface inside the left (right) barrier that separates 
the dot from the corresponding reservoir. It is impor- 
tant to point out that 4>k( x ) m Eq. ([2]) is a bound state 
wave function for the "inner" potential. The latter co- 
incides with the original potential inside the surface E 
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and a constant outside this region. On the other hand, 
Xi( r )( x ) is a non-resonant scattering wave function in the 
"outer" potential, which coincides with the original po- 
tential outside the surface S and a constant inside this 



regionii. 



Taking n along a coordinate of the tunneling electron, 
xn + i, we can integrate over x\, . . . ,xn in Eq. @ thus 
reducing this equation to Eq. ^ with 4> n being replaced 
by the overlap function 



In one-dimensional case (Fig.[T|), the separation surface 
£ becomes the separation point, x, inside the barrier, 
Fig. [1] Then Eq. ^ can be rewritten asi 2 - 
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Therefore, the sign of fl^fl^ is determined by the sign 
of ip N (xi)cp N (x r ). 

By applying the mean-field approximation, we can 



write |<I> 



where = yf2m[UD{xi,r) — Ek] an d (j>k (x) is the bound 
state wave function in the potential Ud{x) (Fig. [1]). The 
separation points xi jr are to be taken inside the left 
(right) barrier as indicated in Fig. [T] and far away from 
the classical turning points^. 

We start with non-interacting electrons, He — in 
Eq. (p}. Then the electron transport through the level 
E\ can be described by the time-dependent Schrodinger 
equation ihdt\^{t)) = H\^>(t)) for a single electron. Tak- 
ing the stationary limit we obtain the Landauer formula 
for the total current, with the transmission amplitude 
given by the Bright- Wigner formula 
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and as a product of one-electron 
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where TV = -2it(qlQm) 1/2 and T^ R = 2^^) V,fl 
are the partial widths, and Ql{R) is the density of states in 



L l,r 



(A) 
L,R 



the left (right) reservoir. We assumed that Q) 
are weakly dependent on E^ r . 

The corresponding evolution of the resonance trans- 
mission phase for different states |A) is determined by 
sign of the product of Q^Q R K Since the reservoir states 
Xi, r are not affected by the plunger voltage, one finds 
from Eq. ^ that the evolution of the sign [fi^ fi^] is 
given by the sign of the product 4>x(xi)(px(x r ). The lat- 
ter is positive or negative, depending on the number of 
nodes of <fi\(x). Hence, it is clear that the non- interacting 
electron model cannot explain the same sign for all reso- 
nances, observed in Refi^ (see also Refs^^). 

Consider N interacting electrons trapped inside the 
dot. Despite the electron-electron interaction, the cou- 
pling amplitudes £Il,r can still be evaluated by us- 
ing the same multi-dimensional overlapping formula 
([2]), as in the non-interacting case. Indeed, the 
many-body tunneling can be considered as one-body 
tunneling, but in the many-dimensional space. In 
this case, the wave-function xi(r)( x ) 1S replaced by 

Xi(r)(xN+i)$N\ x ii ■ ■ ■ ,%n), where Xl{r) is the wave 

function of tunneling electron in the left (right) reser- of the lowest energy state, (xi) = <f> (xi), in the in 
voir and is the ground state wave function of N ner P ot ential U D , Fig.ffl This wave function is nodeless. 
electrons inside the dot. The wave-function 4>k(x) cor- 
responds to $$+i( x i> ■ • ■ ! x n+i), which is the lowest en- 
ergy state (ground state) of N + 1 electrons in the inner 
potential of the dot (Ud in Fig. |T|) . 



states (orbitals) in the effective single-particle potential, 
Ud + Uc, where Ud is the inner part of quantum-dot 
potential (Fig. [T]) and Uc(x) is the mean-field describ- 
ing the electron-electron interaction. As a result, the 
overlap function Lp^(x) is a bound state wave function 
in the potential Ud{x) + Uc(x), corresponding to one 
of the orbitals. Since the lowest energy state is always 
nodeless^, one might assume that <£n{x) is also a node- 
less one, so that the sign of LpN(xi)ipN(x r ) would be the 
same sign for all resonances. This, however, is not correct 
because of the Pauli principle. Indeed, due to the anti- 
symmetrization, any two orbitals in the product of the 
wave functions representing I'l'^x) cannot be the same. 
Since the lowest state is already occupied, the wave func- 
tion y)jv(^Af+i) must correspond to a higher non-occupied 
orbital, and therefore it cannot be nodeless. Hence, the 
Pauli principle would create serious problems in any at- 
tempt to explain the same sign for all resonances 2 in a 
framework of the mean-field description of the electron- 
electron interaction. 

Note that this problem cannot be resolved even by as- 
suming large coupling with reservoirs, so that the reso- 
nances are overlap. Indeed, the problem is related only 
to the inner component of the resonant state, Eqs. @ 
and . The latter is eventually brought by the plunger 
below the Fermi level, /i^, blocking an appearance of 
the resonance above the Fermi level with a similar inner 
component. 

The same situation holds in a more general case, when 
the interaction term Uc varies with each new electron 
trapped inside the dot, Uc — > Uq (Koopman's theo- 
rem is violated). One finds that due to the central sym- 
metry of the self-consistent potential such a variation of 
Uc with N would not affect the number of nodes in the 
overlap function lpn(xn+i)- As a result, the sign of the 
transmission amplitude would fluctuate between ±1 for 
different resonances. 

We illustrate this point by evaluating the overlap func- 
tion (p N (x N+ i), Eq. ([5]), for N — and N = 1. In 
the first case, fo(xi) coincides with the wave function 

(0), 
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The second overlap function is fifa) = (x2,$\ |3> 
where <& 2 ®\x ll x 2 ) = [4>o(xi)4>i(x 2 ) - cf>o(x 2 )cj) 1 (xi)]/^/2 
is the lowest energy state of two electrons in the potential 

_ (2) 

Ud + U c . Here 4>o.i represent the two first orbitals in 
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this potential. One easily finds that 

<Pi(x 2 ) = J <Z> ( i ) (x 1 )<P i °\x 1 ,x 2 )dx 1 = c o 0i(x 2 ) , (6) 

where Co = J 4>Q{xi)4>o{xi)dxi/^/2. (The second term is 
zero, since (fio and <\>\ are orthogonal due to the opposite 
parities). Therefore, ipi contains one node, so that the 
corresponding transition amplitude changes its sign. 

The same behavior of the overlap function would per- 
sist for any TV. For instance, one easily obtains for TV = 2 
that (^2(^3) oc 00^3(2:3)— ci3(pi(x3), where the coefficients 
c o = ((/>i\<Pi), C13 = (0o 102) and 4>, <j> are the orbitals in 
the potentials, Ud{x) + U(P (x) and Ud{x) + Uq (x), re- 
spectively. Since C13 <§C cq, the overlap function if2 would 
contain an additional node in a comparison to tp\. Thus, 
by assuming the ./V dependence of the mean-field effec- 
tive potential, we are still not able to explain the puzzling 
behavior of the transmission phase. 

The above consideration was based on symmetry argu- 
ments applied to electrons moving in a spherically sym- 
metric mean-field central potential. In fact, the central 
mean-field picture for two-dimensional quantum dots was 
challenged in recent publications*'*'*'*. It was suggested 
that due to the strong inter-electron repulsion inside the 
dot, spontaneous symmetry breaking takes place lead- 
ing to the formation of electron molecules. As a re- 
sult, the electrons appear on the ring (rings) around the 
dot's center. This idea was substantiated by unrestricted 
Hartree-Fock calculations or by using other computa- 
tional techniques****!*. 

In principle, if the symmetry is broken, the overlap 
function ipn(xn+i), Eq. ([5]), could be very different from 
the corresponding orbital (j>N(xif+i) m the spherical sym- 
metric potential. Therefore, it is desirable to investigate 
the evolution of the transmission phase in this case. Con- 
sider again the overlap function ipi{x%) — (x 2 , 
but now without the mean-field approximation, as in 
Eq. ([5]). In fact, by taking the parabolic confining poten- 
tial, the two-electron wave function |$ 2 °^) can be exactly 
calculated**, since relative and center-of-mass coordi- 
nates of two electrons are decoupled in the total Hamilto- 

nian. As a result, $2 ( x i,x 2 ) = 4>cm{xx -\-x 2 )4> r (x 2 — x\), 
where 4> r (—x) = —<p r (x) due to the Pauli principle. Such 
a wave function peaks for x\ = — x 2 and therefore it 
would bear the features of a two-electron molecule*. One 
finds from Eq. ([5]), 

Pi{x 2 ) = J 4>a{xi)4> cm {xi + x 2 )4> r {x 2 - xi)dxi . (7) 

Taking into account that the values of x\ which mainly 
contribute to the integral ([7]) are localized inside the dot 
and that the wave function (j) r (x) is the odd one, we find 
that the overlap function changes its sign when the ar- 
gument varies from xi to x r , Fig. [T] Hence, ipi displays 
one node, as in the spherically symmetric mean-field po- 
tential, Eq. ©. 



One can continue with the same arguments for the 
three and more electron molecules, where the electrons 
are placed on the ring. The corresponding ground 
state wave functions ) would represent a fully 
anti-symmetrized product of the original (site) nodeless 
orbitals^*. Yet, the overlap function Eq. © cannot be 
nodeless. As a result, the sign of ipN(xi)<fN(xr) would 
fluctuate with TV. One can demonstrate it rather easily 
for TV = 3,4. Although it would be hard to extend such 
a demonstration for large TV, there is no reason to expect 
that the sign of tp N {xi)^N{x r ) ceases to fluctuate when 
TV increases. 

It seems from the above arguments that the rota- 
tional symmetry breaking (the electron-molecule forma- 
tion) cannot explain the evolution of the transmission 
phase observed in the experiments*!*. Nevertheless, there 
is an additional feature of the electron molecule, which 
has not been yet utilized. That is due to the electrons 
located on the ring (rings) would develop an additional 
(inner) electrostatic trap inside the dot when their num- 
ber (TV) is large enough. As an example, we display in 
Fig. [2k. such a potential, Vc(x) = J^f e 2 /e\x — Xj\, pro- 
duced by 14 electrons equally distributed on the ring, 
where e = 13.6 is the dielectric constant of the medium. 
The radius of the ring (i? — 50 nm) is taken close to the 
dot's size in RefA The radial profile of this potential 
along the angle 9 — tt/N, where the potential height on 
the ring is minimal, is shown in Fig. [2b for two values 
of TV. It appears that the trap is not well developed for 
TV = 6, but it is already pronounced for TV = 14. 

A minimum number of electrons in the dot sufficient 
to develop the trap with one bound state inside it can 
be estimated from the condition that the barrier height, 
hjy in the Fig. [2b, reaches the ground state energy eq. 
We estimate the latter as ir 2 h 2 /m*R 2 , where m* is the 
effective electron mass (to* /too = 0.067). For instance, 
one finds £0 = 4.5 meV for R — 50 nm. Then the condi- 
tion fix — £q corresponds to TV ~ 10, which is a minimal 
("critical") number of electrons, TV cr , enabled to hold a 
resonance state. This value is an approximate agreement 
with that found in*. In fact, a more elaborate, semi- 
classical estimations of N cr approximately produce the 
same number [h N ~ 10 meV for TV = 14, Fig. Wp^- 

The state |e ) in the inner part of the trap, Vc, Fig. [2b, 
is not stable due to the symmetry breaking, leading to 
formation of the electron molecule. Nevertheless, this 
state is important in formation of the (TV + l)-electron 
molecule by adding an additional electron to the TV- 
electron system. Indeed, one expects that the overlap 
function for the electron states on the ring is sup- 
pressed in comparison to the same overlap for the cen- 
tral mean-field potential. The reason is that all electrons 
are shifted from their positions whenever an additional 
electron is placed on the ring. This is in contrast to the 
mean-field description, where the TV-electron core is not 
modified. On the other hand, if the electron is placed in 
the center of the ring, it distorts the remaining TV elec- 
trons in a minimal way. We expect therefore that the 
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FIG. 2: (Color online) (a) Electrostatic trap generated by 14 
electrons placed on the ring of the radius 50 nm. (b) The 
radial profile of the Coulomb potential along the potential 
valley, £o is the ground state energy in the potential Vc, 
representing the inner part of the trap Vc- 



corresponding overlap function is large, as in the case of 
the central mean-field potential. Hence, such an unstable 
state |eo) m the middle of the dot would play a role of 
a "doorway" state in formation of the (N + l)-electron 
molecule. 

It follows from the same arguments that the electron 
transport would proceed through such an unstable state 
when the quantum dot coupled with the reservoirs. Since 
this doorway state is of the lowest energy in the inner 
trap, Vc (Fig- [2]), it is nodeless. The crucial point here 
is that this state is eventually not occupied, when it is 
brought by the plunger below the Fermi levels of the 
reservoirs. Indeed, it is not turned to a stable state be- 
low the Fermi levels due to the symmetry breaking, but 
it always decays to the ring states. Therefore, this state 
is never blocked by the Pauli principle to carry the res- 
onant transport through it, when it is above the Fermi 
level iir, Fig. [I] As a result, all transmission amplitudes 
for any N > N cr would be in phase. 

In fact, by taking into account the electron spin, one 
finds that two electrons with the same spatial (nodeless) 
wave functions are allowed to occupy the lowest energy 
states. Therefore, even if the state \eo) in the center of 
the dot becomes a stable one for some values of N, the 
resonant transport would proceed through an unstable 



state of the two electrons (with opposite spin) inside the 
dot. The corresponding overlap function would be again 
nodeless. 

Ncr 
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FIG. 3: Dependence of the "critical" number of electrons on 
the dot's radius. 



Note that although the doorway-state energy is the 
lowest one for the inner trap, Vc, it exceeds the energy 
of the ring states. Therefore, the ring states would ap- 
pear inside the bias voltage before the doorway state. We 
can assume, however, that the ring states are not well 
separated in the energy from the doorway state, which 
dominates the resonant current. It was also taken into 
account that in the presence of the Coulomb interaction, 
the shift of the resonance energy due to tunneling is dif- 
ferent for different levels^. In particular, the broad reso- 
nance is shifted down more than the narrow onei^. As a 
result the doorway state could have a lower energy than 
the ring states. 

One of the consequences of our model is an existence of 
the critical number of electrons in the dot, which is nec- 
essary for formation of the resonant state inside the dot 
(N cr ). This number would vary with the dot's size. Such 
a dependence of N cr on the radius of the dot (R), ob- 
tained from our estimation, /ijv = £oi is shown in Fig. [3] 
One finds from this figure that this dependence is rather 
weak. The critical number slightly decreases with an in- 
crease of the dot's size. 

In summary, we demonstrated that the unusual be- 
havior of the resonant phase, observed in interference ex- 
periments, can be considered as a strong evidence for 
formation of electron molecules in quantum dots. This 
structure would produce an electrostatic trap, contain- 
ing an unstable (doorway) state localized in the center of 
the dot, whenever the number of electrons occupying the 
dot is large enough, N > N cr . Then such an unstable 
state would carry the electron transport through the dot 
irrespective of the value of N. This would appear as if 
the different transmission amplitudes are in phase. Our 
prediction for the dependence of N cr on the dot's radius 
can be experimentally verified. 
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